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The state matrix p for an open quantum system with Markovian evolution obeys a master 
equation. The master equation evolution can be unraveled into stochastic nonlinear trajectories for 
a pure state P, such that on average P reproduces p. Here we give for the first time a complete 
parameterization of all diffusive unravelings (in which P evolves continuously but non-differentiably 
in time). We give an explicit measurement theory interpretation for these quantum trajectories, in 
terms of monitoring the system's environment. We also introduce new classes of diffusive unravelings 
that are invariant under the linear operator transformations under which the master equation is 
invariant. We illustrate these invariant unravelings by numerical simulations. Finally, we discuss 
generalized gauge transformations as a method of connecting apparently disparate descriptions of 
the same trajectories by stochastic Schrodinger equations, and their invariance properties. 



I. INTRODUCTION 

It is well known that quantum mixtures differ qualita- 
tively from classical mixtures. Here a mixed state means 
one about which we have incomplete knowledge, as op- 
posed to a pure state, which is one about which we have 
maximal knowledge. In the classical mixed state 

is described by a probability distribution p over phase 
space. There is a one-to-one correspondence between 
this probability distribution p and a weighted ensemble 
of points in phase space (pure states). In the quantum 
case, a mixed state is described by a density operator, 
or state matrix p. But there is a one-to-many (infinitely 
many, in fact) mapping p) from this p to a weighted 
ensemble of state vectors (pure states). That is to say, 
there are infinitely many ways to write a given impure 
state matrix p as a positively weighted sum of projectors. 

This difference between quantum and classical systems 
is also reflected in the dynamics of open systems. Inter- 
action with an environment generically causes systems to 
become mixed. This can be described by deterministic 
evolution of the mixed state (the classical p or the quan- 
tum p). For example, this evolution may be described 
by a Fokker-Planck equation for p, or a master equation 
for p. Alternatively, the dynamics can be described by 
stochastic trajectories of pure states (the classical point 
in phase space or the quantum projector). In the classical 
case the two descriptions follow uniquely from each other 
and they are completely equivalent. But in the quantum 
case, the trajectory equation does not follow uniquely 
from the ensemble equation. In fact there are infinitely 
many different quantum trajectory equations for a given 
master equation |^||. 

It is thus apparent that quantum trajectories have a 
richer physical content than classical trajectories. The 
different quantum trajectory equations have been called 
different unravelings of the master equation. In some 
(but not all Q]) cases, they can be interpreted as arising 
from inequivalent schemes of efficiently monitoring the 
environment to which the system is coupled || . (For clas- 



sical systems, all such efficient schemes would be equiva- 
lent with one another.) From this point of view, the sta- 
tistical description arises from the probabilistic nature of 
quantum measurements. However, it is possible to study 
the nature of quantum unravelings without specifying a 
concrete measurement scheme. 

In the present work, we limit our discussion to a par- 
ticular class of unravelings, in which the noise in the 
quantum trajectory is diffusive in nature. After some 
preliminaries, we give, for the first time, a complete pa- 
rameterization of such unravelings in Sec. IV. An explicit 
formulation in terms of general quantum measurement 
is given in Sec. V. In Sec. VI we discuss the notion of 
invariance of quantum trajectories, using concepts from 
Sec. II. We introduce new classes of invariant unravelings 
and, in Sec. VII, illustrate them for the process of reso- 
nance fluorescence. In Sec. VIII we discuss the stochastic 
Schrodinger equation formulation of quantum trajecto- 
ries, and show how seemingly inequivalent equations may 
be related by generalized gauge transformations. We con- 
clude with a summary and discussion of the relation of 
our present work to past and future work in the field. 



II. THE MASTER EQUATION 

If a quantum system is weakly coupled to an environ- 
mental reservoir, and many modes of the reservoir are 
roughly equally affected by the system, then one can 
make the Born and Markov approximations in describing 
the effect of the environment on the system || . Tracing 
over (that is, ignoring) the state of the environment leads 
to a Markovian evolution equation for the state matrix 
p of the system, known as a quantum master equation. 
The most general form of the quantum master equation 
which is mathematically valid is the Lindblad form j?J 

P = Cp = -i[H,p\ +c k pc\ - | j4c fc ,pj . (2.1) 



f 



Here {ck}^ =1 is the ordered set of Lindblad operators, 
and as in the remainder of this paper, we are using the 
Einstein summation convention for repeated indices. The 
Lindblad operators couple the system to the reservoir 
modes. 

The above representation of the evolution superoper- 
ator C is not unique. We can reduce the ambiguity by 
requiring that the operators 1, c±, c%, . . . , ck be linearly 
independent. Then we are left with the freedom of re- 
defining the Lindblad operators by an arbitrary K x K 
unitary matrix T^i H : 



Cfc 



Tkici . 



(2.2) 



Here, as in the remainder of this paper, we are using 
the Einstein summation convention for repeated indices. 
In addition, C is invariant under c-number shifts of the 
Lindblad operators, accompanied by a new term in the 
Hamiltonian 



Cfc -> c fc + Xfc 



H 



ff+i(XfcCfc-H.c.) . (2.3) 



The master equation turns pure states p = \ip) (ip\ = P 
into mixed ones. A related mathematical object is the 
transition (mixing) rate operator M 



W = CP- {P, CP} + PTi [PCP] 

= (Cfc - (cfc))P(Cfc- (Cfc)) f , 



(2.4) 
(2.5) 



where (...) ~ Tr[. . . P] stands for quantum expectation 
values. W is an invariant op erator; it does not change 
with the unitary rotati on (|2.2| ) of the Lindblad operators 
or the c-numbcr shift ( |2.3p . The transition rate opera- 
tor is positive semi-definite and orthogonal to the current 
state P, i.e. WP = PW = 0. Its trace is the transition 
(mixing) rate: 



TrW 
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fc C fc/ 



Cfc, 



(2.6) 



III. QUANTUM TRAJECTORIES 

In the situation where a Markovian master equation 
can be derived, it is possible (in principle) to continually 
measure the state of the environment on a time scale large 
compared to the reservoir correlation time but small com- 
pared to the response time of the system. This effectively 
continuous measurement is what we will call "monitor- 
ing" . In such systems, monitoring the environment does 
not disrupt the system-reservoir coupling and the system 
will continue to evolve according to the master equation 
if one ignores the results of the monitoring. 

By contrast, if one does take note of the results of mon- 
itoring the environment, then the system will no longer 
obey the master equation. Because the system-reservoir 
coupling causes the reservoir to become entangled with 
the system, measuring the former's state produces in- 
formation about the latter's state. That is, the system 



state is conditioned upon the result of the measurements. 
This will tend to undo the increase in the mixedness of 
the system's state caused by the coupling. 

Perfect monitoring of the reservoir requires continual 
rank-one projective (i.e. von Neumann) measurements 
of its state, on the time scale discussed above. If the sys- 
tem initially has a mixed state, then its state will usually 
be collapsed towards a pure state. However this is not a 
process which itself can be described by projective mea- 
surements on the system, because the system is not being 
directly measured. Rather, the monitoring of the envi- 
ronment leads to a gradual (on average) decrease in the 
system's entropy. 

If the system is initially in a pure state then, under 
perfect monitoring of its environment, it will remain in 
a pure state. Then the effect of the monitoring is to 
cause the system to change its pure state in a stochas- 
tic and (in general) nonlinear way. Such evolution has 
been called a quantum trajectory p0| . It can be de- 
scribed by a nonlinear stochastic Schrodinger equation 
(SSE) @||e,|,p|-P]: 



= —iHjf, \ip) + noise,/, 



(3.1) 



Here is a non-Hermitian effective Hamiltonian, and 
(like the noise) it depends on ip. The nonlinearity and 
stochasticity are present because they are a fundamental 
part of measurement in quantum mechanics. 

The stochastic average of pure state quantum trajec- 
tories still reproduces the state of the ensemble p for each 
time t: 



e[ m))W)\ \=E[p(t)] = p{t) 



(3.2) 



Here E denotes an expectation value, or ensemble aver- 
age with respect to the noise process in the stochastic 
Schrodinger equation. The non linear Hamiltonian and 
the stochastic term in Eq. ( |3.l| ) must be derived from 
the above constraint. Then the stochastic Schrodinger 
equation is said to unravel the master equation It 
is now well-known |l6| ] that there are many (in fact con- 
tinuously many) different unravelings for a given master 
equation, corresponding to different ways of monitoring 
the environment. 

Any classical noisy trajectory can be approximated to 
an arbitrary accuracy by a trajectory consisting of de- 
terministic evolution punctuated by jumps |JLfl|. In the 
same wa y, the noise in the stochastic Schrodinger equa- 
tion (3^) can always be written as a quantum jump term. 
These jumps may range in size from being infinitesimal, 
to being so large that the system state after the jump 
is always orthogonal to that before the jump P, |l8| , |l9| . 
In this paper we are concerned with the former case. In 
the limit of infinitesimal jumps occurring infinitely fre- 
quently, a diffusive unraveling results. As in Brownian 
motion, the state of the system evolves continuously but 
not differentiably in time. For this reason, these sorts of 
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unravelings have been called continuous unravelings pl| , 
but here we call them diffusive. 

Although a stochastic Schrodinger equation is concep- 
tually the simplest way to define a quantum trajectory, in 
this work we will instead use the stochastic master equa- 
tion (SME) H§C] This has a number of advantages. 
First, it is more general in that it can describe the purifi- 
cation of an initially mixed state. Second, it is easier to 
see the relation between the quantum trajectories and the 
master equation which the system still obeys on average. 
Third, it is invariant under gauge transformations 



\dcp) (d<p\ = Wdt 



(4.5) 



(3.3) 



where x(f) ^ s an arbitrary real function of time. Such 
gauge transformations can radically change the appear- 
ance of a stochastic Schrodinger equation, since x{t) may 
be stochastic and may be a function of \ip) itself. Since 
unnormalized w ave f unctions may also be used, the gauge 
transformation ( |3.3| ) can be extended for complex func- 
tions We discuss these points further in Sec. VIII. 



IV. DIFFUSIVE UNRAVELINGS 

Assuming that the initial state of the system is pure, 
the quantum trajectory for its projector will be described 
by the SME of the form P = CP + noise. The drift term, 
by virtue of Eq. ( |3.2[), a ssures the consistency with the 
ensemble evolution ( |2.l| ). The noise term, which has an 
expectation value of zero, we are assuming to be diffu- 
sive in nature. It is convenient to represent such singular 
Gaussian noise using the Ito calculus (z| . Writing its Ito 
-differential as dX, we have the following general form of 
SME: 



dP = dtCP + dX 



(4.1) 



Here dX is traceless and Hermitian, and depends nonlin- 
early on the current pure state P. Its expectation value 
E[dA] is zero. Since the SME is assumed to preserve the 
purity of the state, the second moments of dX are con- 
strained by the identity dP = d{P 2 ) = {dP, P} + dPdP, 
which must hold for arbitr ary rank-one projectors P. 
Substituting the expression (4.1) and using the Ito rules 
yields two separate equations: 



dX = {P, dX} , 



dXdX = CPdt - {P, CP} dt 



(4.2) 



(4.3) 



We can write the general solution of Eq. ( |4.2| ) in the 
simple form: 



dx = \dv) (ip\ + 1</>> (cM , (dip\ tP) = o 



(4.4) 



Here \d<p) is an Ito -differential of zero mean, orthogo- 
nal to the current state liM- Its autocorrelation is con- 



where W is the P-dependent transition rate operator 
( |2.4|) . Hence we have obtained the general form of diffu- 
sive unravelings in terms of the following SME: 



dP = dtCP + \dip) + \tp) (dip\ , 



(4.6) 



where the Ito -differential \d<p) is orthogonal to \ijj), has 
zero mean, and the Hermitian part of its correlation is 
given by Eq. (4.5). 

The careful reader will observe that the non- Hermitian 
part \dtp) \dip) remains free, expressing the fact that there 
are infinitely many pure state diff usiv e unravelings of the 
same Lindblad master equation (2.1). If the correlation 
\dip) \dip) is set to zero then the noise term is uniquely 
defined by the Hermitian correlation in Eq. (^^), and we 
obtain a unique unraveling 1 23 24f] . Let us call it the stan- 
dard one. The standard SME ( [4.6|) follows uniquely from 
the Lindblad master equation Q2.1| ) and it is invariant in 
a sense that it does not change with the re-definition 
(pE2 , 2.3) of the Lindblad operators. For a long time it 
has apparently been thought that the standard one is 
the only invariant unraveling p3|,ppS[]. In Sec. VI, how- 
ever, we display invariant and non-zero choices for the 
non-Hermitian correlation \dip) \dtp) of the noise. 

Equation (4.6), with the constraints listed below it, 
represents diffusive unraveling in complete generality. 
However, for many practical purposes, it is useful to 
have a more explicit construction. That is we wish to 
reparametrize the state- vector valued generalized Wiener 
noise \dtp) by complex- number- val ued standard Wiener 
noises. Recall the representation ( |2.4| ) of the transition 
rate operat or W in terms of the Lindblad operators c k . 
From Eq. (4.5) it is then obvious that \dip) is spanned 
by the vectors (ck — {c k )) We introduce the vector 



of complex Wiener processes £(t) = {£&(*)} 
cients: 

\dip) = ( Ck -{ Ck ))\i>)de k . 



K 
k=l 



as coeffi- 



(4.7) 



Recall that we are using the Einstein summation conven- 
tion, and note our notation of using £(f) for the Wiener 
process, not its time-derivative as in Ref. p2] j. 

The mean increments E [d^ k ] vanish, so that E [\dip)] 
does als o. A lso, the above form of \df) satisfies the con- 
straint (4.5) provided the Hermitian part of the noises' 
correlation matrix is the unit matrix, while the non- 
Hermitian part remains free: 



d^(t)dC k (t) = dtd Jk 

d£j(t)d£ k (t) = dtuj k 



(4.8) 
(4.9) 



strained: substituting Eq. (4.4) into Eq. (4.3) yields 



The Uj k — u k j are arbitrary complex numbers subject 
only to the condition that the 2K x 2K correlation ma- 
trix of the real vector ^Rec?^*, Imd^j 
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dt 
T 



1 + Re[u] 
Im[u] 



Im[u] 
- Re[u] 



(4.10) 



be nonnegative. It can be shown that the smallest eigen- 
value of this matrix is given in terms of the norm of the 
complex matrix u by (1 — ||u||)/2. Thus, this condition 
is satisfied if an only if 



(4.11) 



lull < 1 



With this parameterization we can rewrite the SME (|4 
explicitly as 



dP = CPdt + [( Cfc - (c fc » PdC k + H.c] 



(4.12) 



V. MEASUREMENT INTERPRETATION 



We stated in Sec. Ill that the master equation is un- 
raveled if the environment of the syste m is monitored, 
and that the pure state P obeying Eq. ( 4.12 ) can be in- 
terpreted as the state conditioned on the results of this 
monitoring. To see this relationship, it is necessary to 
consider the theory of non-projective or indirect mea- 
surements (see for example Ref. pEj] ). Such measure- 
ments arise when the system of interest interacts with a 
second system, and that second system is subject to a 
measurement of the traditional (projective) sort. If the 
second system is initially in a pure state, and a rank-one 
projective measurement is made on it, then the indirect 
measurement on the system can be described by a set 
of measurement operators Q r . Here r labels the result of 
the measurement, and the operators are constrained only 
by the completeness relation 



1 . 



(5.1) 



Here dpo (r) is a normalized measure over the space of all 
r. Let the initial state of the system be p. The measure- 
ment operators give both the probability 



dp(r) — d^o( r )Tr [pO^Or] 



(5.2) 



for obtaining a result in an infinitesimal vicinity of r, and 
the state 



dfio(r)n r p£ll tt r pQJ, 



dn(r) 



Tr 



pfttflj 



(5.3) 



conditioned on the result r. If the measurement is made 
but the result ignored, the new system state is 



p' = dfi(r)p' r = / dfio(r)Sl r pdl 



(5.4) 



In this paper we ar e con cerned with systems that obey 
the master equation (2.1), which is continuous in time. 



That is to say, we have to consider repeated indirect mea- 
surements, each lasting an infinitesimal (with respect to 
the relevant system time scales) interval of time, such 
that if one ignores the result, one obtains 



p = p + dp = p + dtCp . 



(5.5) 



In order to obtain the conditioned evolution equation de- 
scribed by the SME (4.12) it turns out that the measure- 
ment result r in any infinitesimal time interval [t, t + dt) 
must be described by a vector of complex numbers J(t) = 

{Jk(t)}k=i 4 As functions of time, these are continuous 
but not differentiable, and we will call them currents. 
Explicitly, the y are related to the complex Wiener incre- 
ments in Eq. (4.12) by 



Jkdt = dt ( UkjC 



Ck 



d£k 



(5.6) 



That is, it is the randomness in the measurement record 
which provides the stochasticity in the quantum trajec- 
tory. 

We can prove this relation between the noise in the 
quantum trajectory and the noise in the measurement 
record by using the theory of indirect measurements de- 
scribed above. We define the measurement operator to 
be 



iHdt 



k Ckdt + Jlckdt 



These obey the completeness relation 



d/i (J)OtOj-= 1 
if we choose dpo(J) to be the measure such that 
dp (J)(J k dt) = 



(5.7) 



(5.8) 



o , 


(5.9) 


Sjkdt , 


(5.10) 


Ujkdt . 


(5.11) 



dp (J){J*dt){J k dt) 



dpo(J)(Jjdt)(Jkdt) = Ujkdt 



These moments are the same as those of the Wiener in- 
crement dt; as defined above. 

With this assignment of measurement operators £1 and 
measure dpo we can easily show that the expected value 
of the result J is 



E[J fc 



dp{J)J k = lukjc] + c fe 



(5.12) 



This is consistent with the previous definition in 
Eq. (5.6). Furthermore, we can show that the second 
moments of Jdt are (to leading order in dt), indepen- 
dent of the system state and can be calculated using d/io 
rather than dp. In other words, they are identical to 
the statistics of d£ as defined above. This completes the 
proof that Eq. ( |5.6D gives the correct probability for the 
result J. 
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The next step is to derive the conditioned state of the 
system after the measurement. According to the theory 
of indirect measurements this is given by 



and 



P + dP = 



Tr 



fljPflL 



(5.13) 



Expanding to order dt gives 
dP = -\ |4c fc) P} dt + Jldtc k P ciJidt 
+{JtdtJk - 1) (c\cj\ Pdt 

+J* k {ck - (c k ))Pdt + P{c\ - {ck)*)J k dt . (5.14) 



Substituting in t he ab ove result fl5.q ) for J yields the re- 
quired equation (4.12). From this it is again obvious that 
on average the system obeys the master equation. In the 
measurement interpretation this can be derived directly 
from the nonselective (ignoring the measurement result) 
evolution 



dp = J dn{J)dP = -P + J dno(j)Cl f Pnt . (5.15) 

Some insight into the above formalism may be found 
by considering an experimentally realisable situation in 
quantum optics [p6f . For simplicity we consider a system 
with one irreversible term; that is, K = 1. For speci- 
ficity, say the system is a two-level atom, with sponta- 
neous emission rate 7. Then c = y/jcr, where er = \g) (e| 
is the lowering operator for the cavity. Say the atom is 
placed at the focus of a parabolic mirror so that the emit- 
ted light emerges as a beam, and let that beam impinge 
upon a beam splitter of transmittance 77. Let the trans- 
mitted beam be subject to homodyne detection with a 
local oscillator of phase (relative to the system) of 0\. 
This will yield a real homodyne photocurrent of |Tc| ] 



hdt 



jrj(e- iei a + e idi ^) + dCi{t) , 



(5.16) 



which has been normalized to have a shot-noise spectrum 
of unity from the real Wiener process (i(t). Let the re- 
flected beam be subject to homodyne detection with a 
local oscillator of phase #2, yielding 



I 2 dt = v/ 7 (l - 77) (e^a + e"***) + d£,(t) , (5.17) 

where is an independent real Wiener process. 

From these real currents we can define a 'complex cur- 
rent' 



Jdt = y^ie i9l hdt + y/l - r]e i02 I 2 dt 
= (u^ + °) dt + d((t) , 



where 



2,1)-, 



(5.18) 
(5.19) 



(5.20) 



d£ = ^oe^dd + ^n^e^db 



(5.21) 



obeys d£*d£, = dt, d£d£ = udt. Furthermore, the condi- 
tioned system state can be shown to obey the expected 
stochastic master equation ]27|. 



VI. INVARIANT DIFFUSIVE UNRAVELINGS 



From Eq. (4.12) it is obvious that all diffusive unrav- 
elings with u fixed are invariant under the shift transfor- 
mation of Eq. (2.3). However, from Eq. ( |4.9D we see that 
the unraveling is in general not invariant against the uni- 
tary rearrangement (2.2) of the Lindblad operators. It 
is invariant if the non-Hermitian correlations ujk vanish 
(the standard diffusive unraveling). As noted above, it 
has been guessed that this was the only invariant un- 
raveling. We show here that there are further invariant 
unravelings. 

Consider the following non-Hermitian correlations 



(4.9) 



u jk = R x ((cj - (cj)) (c fe - (c fe ))) 



(6.1) 



Here R is a complex number constrained only by the fact 
that its magnitude must be sufficiently small for the pos- 
itivity condition (4.11) to be satisfied (note that ||u|| is 
invariant). For a system with unbounded Lindblad oper- 
ators c/c, the invariant number R would have to depend 
upon P to ensure this. An obvious choice would be for R 
to be real, and to take the maximally po sitive (or max- 
imally negative) value that satisfies Eq. ( 4.11 ). For the 
special case of finite A^-dimensional Lindblad operators, 
state independent alternatives can also be chosen, for ex- 
ample 

u jk = RxTt [(cj - N~ lr TrCj) (c k - AT 1 !*^)] . (6.2) 

The abo ve correlations are trivially invariant for the 
shifts (2.3). It is crucial to notice that the coefficient Ujk 
depends on operator product CjCk instead of the Hermi- 
tian versions Cjc\ or c^Cfc. This little difference assures 



that the SME (4.12) will be invariant for rotations (gj). 



To inspect this invariance, observe that the Lindblad op- 
erators and the co mple x noises appear always in the same 
combination \d<p) ( |4.7| ). The mathematical characteriza- 
tion of this Ito-d iffcrcntial is fully given by the Hermitian 
correlation ( |4.5| ) which is invariant (since W invariant it- 
self) and by the non-Hermitian correlation 

\dtp) ® \dtp) = u* k (cj - (cj)) <g> (c fe - (cfc)) dt . 

(6.3) 



This latter becomes invariant for unitary rotations ( |2.2| ) 
if we use the non-trivial definitions (6.1) or ( |6.2| ) for Ujk- 
It is interesting to note that the choice of correlations 
(|6.l[) implies that the noise process d£(t) is no longer 
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white. That is because the noise correlations depend 
upon the state of the system at that time, which ob- 
viously depends upon past values of the noise. Never- 



theless, the quantum trajectory defined by Eq. (3.1) is 
still Markovian, in that dP depends only upon P at the 
present time, and the noise process d£ is still uncorrelated 
with P. There are other choices of u which would make 
the quantum trajectory strictly non-Markovian. For ex- 
ample, u could depend upon the past values of the cur- 
rent J, and in fact there are very practical reasons for 
wishing to consider such unravelings p8[ . However we 
will not be concerned with this possibility here. 



steady state of the master equation is close to a com- 
pletely mixed state, so that x (exactly), and y and z 
(approximately) average to zero. 

We begin with two non-invariant unravelings, u — 1 
and u = — 1. These correspond to homodyne detec tion 
as described in Sec. |v|. The noise correlations ( 4.10 ) de- 
generate and a standard real white noise C will drive the 
quantum trajectories. For u = 1, the current ( |5.6| ) be- 
comes real, with average 



E[J(*)] = Vl (<r x 



(7.4) 



VII. SIMULATION OF UNRAVELINGS 



and noise £ 
follows: 



(. The SME fl4.12j ) can be written as 



In this section we illustrate various unravelings (invari- 
ant and non-invariant) for a simple but interesting quan- 
tum optical system: a driven, damped two-level atom. 
The master equation in the interaction picture is Iq] 



P = 1°P^ - \ {^e,p} - i4[vx,p] 



(7.1) 



Here fi is the driving strength (the Rabi frequency) and 
damping occurs through spontaneous emission with the 
single Lindblad operator ^Ffa as described in Sec. [v|. 
Physically, all of the light emitted by the atom would 
have to be collected and measured by homodyne-like 
measurements (as described in Sec. ^) in order for a dif- 
fusive SSE to describe the conditioned dynamics of the 
system. 

Since there is only one Lindblad operator, all diffu- 
sive unravelings are defined by just one complex param- 
eter u satisfying |u| < 1. In all cases illustrated below 
the damping rate 7 is set to unity, and the driving rate 
£1 = IO7. We plot the state using the three components 
of the Bloch vector x, y, z, defined by 



P = i (1 + xa x + yo y + za z ) . 



(7.2) 



Here a x ,a y , a z are the usual Pauli pseudospin matrices 
where the up (|e)) and down (\g)) states are the a z eigen- 
states with eigenvalues +1 and —1, respectively. One of 
the Bloch vector components is redundant since for a 
pure state P 2 = P which implies 



V 



1. 



(7.3) 



Nevertheless, it is easiest to interpret the results if we 
plot all three. The driving causes the Bloch vector to 
rotate around the x axis, while the damping causes it 
to decay towards the down state (z = — 1). In all cases 
the initial state is a positive x eigenstate. We emphasize 
that the ensemble average behaviour is the same for all 
unravelings. In particular, after transients have decayed 
the system on average reproduces the stationary solution 
of the master equation. In the high driving limit, the 



dP = CPdt + %[<r y ,P]dt + \{a x - {a x ) , P}d( . (7.5) 

The third term on the SME's RHS corresponds to a mea- 
surement of the x quadrature of the syst em, as reflected 
in the expectation value of the current (7.4). However, 
there is a second white noise term on the RHS, corre- 
sponding to a noisy Hamiltonian. It can be interpreted 
as an additional (non-Heisenberg) back-action caused by 
the monitoring, as if the current noise was being fed-back 
to alter the system dynamics. The presence of two (cor- 
related) noise terms in the SME, one describing Heisen- 
berg back-action and one describing Hamiltonian noise, 
is a generic feature of unravelings with a non-Hermitian 
Lindblad operator (in this case, a). 





1 ff^v 








FIG. 1. Non-invariant unraveling with u = 1 (homodyne 
measurement of x quadrature). 

Fig. shows the conditioned evolution for u = 1. Mon- 
itoring the x quadrature tends to make x well-defined 
(i.e. close to the a x eigenvalues of ±1). However, it is 
certainly not perfect in this respect, as large oscillations 
in y and z due to the rotation around the x axis at rate 
Q, are still evident. 



G 





FIG. 2. Non-invariant unraveling with u = — 1 (homodync 
measurement of y quadrature). 

Figure [2] shows the case u = — 1, which corresponds 
to a homodyne measurement of the y quadrature of the 
system, so 

E[J(t)] = K) ■ (7.6) 

The measurement tries to make y well-defined, but fails 
because the Rabi cycling rotates y into z and so on. Nev- 
ertheless, x is forced towards zero, where it stays. 




FIG. 3. Invariant unraveling with u = (heterodyne de- 
tection). 

The standard invariant unraveling, or "quantum state 
diffusion" || case of u — is shown in Fig.||. This could 
be realized by heterodyne detection pl| , which is like ho- 
modyne detection but with a far-detuned local oscillator. 
This ensures that both quadratures are sampled equally. 
The current J (which is the complex Fourier transform 
of the physical photocurrent) has a mean 

E[J(t)] = ^ <*> • (7-7) 

The resultant evolution is intermediate between that of 
u = 1 and u = — 1: both <j x and a v are being equally 



monitored, and the result is controlling a certain dynam- 
ical feed-back |p9fl . 

Next we plot the simulation of the first of our ne w in - 
variant unravelings. We choose u to be given by Eq. ( |6.l| ), 
where R is chosen to be the p osit ive number such that 
\u\ = 1 [the alternative, Eq. ( |6.2| ), simply gives u = 
again]. In this case we find u = — (a) / (c'), so that 



E[J(t)] = 



(7. 



That is, the current consists purely of white noise. This 
is not because the monitoring no longer gives any infor- 
mation about the system; it is merely because the the 
measured quadrature of a is always chosen to be the one 
that has an instantaneous mean value of zero. 




FIG. 4. Invariant unraveling with u given by Eq. (6.1) with 
maximally positive R. 

The resultant trajectories are shown in Fig. |[ It ap- 
pears that the behaviour is qualitatively similar to that 
of u = 0. However, note that unlike any of the other 
cases, the evolution of z is diffcrentiable. This can be 
proven analytically, since in this case z obeys 



i = — 72 — 7 + Qy 



(7.9) 



which is exact ly th e same as that obtained from the mas- 
ter equation (7T). The evolution of z is of course still 
stochastic because of the coupling to y. This feature is a 
special case of a general phenomenon noted in Ref . Q , 
namely that the noise correlations may be chosen to make 
the evolution of any given operator average smooth. 

The final plot, Fig. |^, is for the new invariant unrav- 
eling with u given by Eq. (3.1), where now R is chosen 
to be the negative number such that |u| = 1. This gives 
u = (a) I (a') so that 



E[J(t)] = V72<<7> 



(7.10) 



In this case the behaviour appears qualitatively similar 
to that of u = —1, in that x is forced to zero. Moreover, 
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the two unravelings becomes equivalent once x reaches 0, 
since then (a y ) = 2(a). This is interesting, in that an 
invariant unraveling and a non-invariant unraveling are 
actually identical in the steady state. 



-0.5 - 

-1 — 





FIG. 5. Invariant unraveling with it given by Eq. (6.1) with 
maximally negative R. 



VIII. STOCHASTIC SCHRODINGER 
EQUATIONS 



We c an r ewrite the general SME (4.6) into the form of 



a SSE (3.1 



d |V) = -iH^dt | V) + (c fc - (c fe )) dCu 



.1) 



Using the identity dP = (d\ip)) (ip\ + (d (ip\) + 
d\ip)d(ip\, one can inspect that the above SSE leads to 
the SME if the nonlinear Hamiltonian is chosen as it fol- 
lows: 



= -iH 



2 (c fc )* c fc 



+ (c fe ) (c fc ; 



(8.2) 
(8.3) 



It is remarkable that H^,, apart from irrelevant c- number 
terms, does not depend on the particular represent atio n 
of the master equation. This is obvious from Eq. (3.2). 
Let us re-group its terms in the following ways: 



H^=H+\(i (c k )*c k +H.c.) 
-| (c fe - (c fc )) f (c fc - (c k )) 



(8.4) 



The second term on the r.h.s. is a nonlinear Hermi- 
tian term, a kind of mean-field potential. It is invariant 
unde r the rotation ( |2.2[ ) and its change under the shift 
( |2.3[ ) exactly cancels the change in H. The third term 
is nonlinear and non-Hcrmitian, it is responsible for the 
localization of the wavefunction as a result of continuous 
monitoring. It is obv ious ly invariant under both trans- 
formations (|2.2|) and (2.2). 



The invariance of is a nice feature of the SSE (|B.l|), 
but it must be emphasized that it is not a required prop- 
erty of a SSE. That is because the mapping from SME 
to SSE is not unique, since a gauge (global phase) trans- 
formation 



(8.5) 



does not change the projector P — \ip) (ip\ ~ \(j>) (<j)\. For 
a non-stochastic Schrodinger equation, this transforma- 
tion will cause only trivial changes to its equation of mo- 
tion; it will simply add a c-number to the Hamiltonian. 
However for a SSE it can radically change the equation. 
This has been noted before (see for example Ref. filf), 
but for completeness we show it explicitly. 

Consider the case K = 1 so that we have one complex 
noise process d£, and one Lindblad operator c which we 
will take to be Hermitian and write as x. Then the SSE 
( |S.1| ) becomes 



d|V>= -iH- \(x- (x)) 2 db\i>) + {x- (x))d£ 



(8.6) 



Now let the global phase obey the equation 

dx = fWC + c.c. , (8.7) 
where f(t) is an arbitrary smooth function of time. Then 

\<j>) + d\<f>) = {l + idx-\d x d X )^ x{t) m + d\^)) . 

(8.8) 

The resultant equation for \<p) is 



d\(j>) = [-iH-Rc(fv 
-| (x- (x)) (x- 



+ \f\ 2 )]dt\t) 

(x) +ifu* + if*)dt\ 



+ {( x -{ x )+if)dt* + if* <%}]<!>) 



(8.9) 



Clearly the deterministic part of this is different from 
Eq. ( |S.6p , and not invariant. 

In this case the transformed equation seems far less ap- 
pealing in form than the original. How ever, one can find 
different forms of the general SSE ( |S.l| ) which, while not 
having an invariant H^, 1 have other attractions. In par- 
ticular, consider the non- unitary gauge transformation 
defined by 

\m) - e ix(t) \m) = \m) , (s.io) 

where \ is the complex function obeying 

id X = (cj) dCj - \u* k (cj) (c fe ) dt+l (cj)* (cj) dt . 

(8.11) 

The state |(/>) is not normalized, which is why it is indi- 
cated with an overbar. This is not important, however, 
as the projector can be defined as P = \cj>) {(j>\ / ((f>\ <f>), 
and this will still obey the SME. 
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Following the same method as above, one finds that 
5) obeys the SSE 



d\<f>) =dt i-iH 



\c\c k + JfeCfc) \(j>) 



.12) 



where it is the normalized state P that is used to define 
the quantum averages in the expression ( |5.6p for the cur- 
rents Jfe. This SSE has a number of nice features. First, 
it is as simple as the invariant version (8.1). Second, it 



very clearly shows the conditioning of the state on the 
measurement result, and is closely related to the mea- 
surement operators (|5.7| ). Third, it is an easy form to 
use for numerical calculation (and was in fact used for 
the simulations in Sec. VII). 



IX. DISCUSSION 

In this paper we have presented new results, and cor- 
rected and clarified old results, pertaining to diffusive 
unravelings of Markovian quantum system dynamics. 

First, we have given for the first time the most general 
form of diffusive quantum unravelings, in Sec. IV. While 
Gisin restricted his early work || for the two-dimensional 
special case, there have been recent publications which 
claim to do much the same thing, but in fact do not do so. 
For example, the recent work of Adler |?2| is restricted to 
Lindblad operators that are Hermitian. Dorsselaer and 
Nienhuis |33| give a general SSE for a master equation 
with one Wiener noise term d£, but for generalizing to 
the set {d£k}k=i they say "we have to assume that the 
different d£k are uncorrelated" . As we have shown here, 
this is not a necessary assumption; d^jd^k — uj^dt may 
be nonzero for j =/= k. The non-Hermitian correlation 
matrix u was introduced by one of us and Vaccaro |n] , 



where it was stated that \fj, k\iijk\ < 1- While this is 
a necessary condition, it is not sufficient. Here we have 
shown that a necessary and sufficient condition is that 
the norm ||u|| be bounded above by unity. 

Second, we have given for the first time the relation 
between the most general unraveling, as parameterized 
by u, and quantum measurement theory. The measure- 
ment results which condition the system and so unravel 
the master equation are K complex "currents" (contin- 
uous functions of time) given explicitly in Sec. V. The 
measurement interpretation of the diffusive unravelings 
is significant in that it means that the ensembles of pure 
states resulting from the unraveling can be physically re- 
alized. Physical realizability was proposed in Ref. pjj ] as 
one of the requirements (along with maximal robustness) 
for finding the "most natural" ensemble of pure states to 
represent the mixed state of an open quantum system. 
Our present work is significant for this program (contin- 
ued in Ref. Q) of investigating decoherence in that gives 
a simple boundary (||u|| < 1) to the parameter space {u} 
of all diffusive unravelings. 

Third, we have corrected the long-standing concep- 
tion p3PJl9|] that the only invariant unraveling is the 



standard diffusive unraveling with u = 0. Here invari- 
ant means invariance under the linear transformations 
of the Lindblad and Hamiltonian operators which leave 
the master equation invariant. In Sec. VI we constructed 
some explicit examples of invariant unravelings with non- 
zero u. The most natural ones have a state-dependent 
u satisfying ||u|| = 1. We illustrated two of these, along 
with the standard invariant unraveling and some non- 
invariant unravelings, by numerical simulations of reso- 
nance fluorescence in Sec. VII. One of the new schemes 
produced quite distinctive dynamics for the atom, which 
ties into recent work on minimizing statistical errors in 
ensemble average simulations using quantum trajectories 
pOf . "Quantum state diffusion theory" g suggests the 
standard unraveling as the "most natural" . The exis- 
tence of non-standard invariant unravelings calls for ad- 
ditional arguments. Recently one of us and Kiefer |35| l 
have app lied robustness criteria (different from those in 
Ref. |2l| ) to open system unravelings and have approved 
the standard one. 

Fourth, in Sec. VIII we have clarified the notion of 
invariance in the context of stochastic Schrodinger equa- 
tions (SSEs) as a way of representing quantum trajecto- 
ries. It turns out that gauge transformations can rad- 
ically alter the structure of a given SSE. In particular, 
the invariance of the standard unraveling is completely 
destroyed in a generic gauge. This suggests that the best 
way conceptually to represent quantum trajectories is as 
a stochastic master equation (SME) for the state projec- 
tor rather than a SSE for the state vector. Gauge free- 
dom may, on the other hand, allow for equivalent SSEs 
with different attractions, as demonstrated. The SSEs 
may be given priority over the SME in numerical calcu- 
lations, but it must be ensured that all theoretical claims 
do not rely on non-gauge-invariant properties. 
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